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Abstract 

The turbulent energy flow of the one-dimensional Maj da-McLaughlin- Tabak equa- 
tion is studied numerically. The system exhibits weak turbulence for weak driving 
forces, while weak turbulence coexists with strongly nonlinear intermittent collapses 
when the system is strongly driven. These two types of dynamics can be distin- 
guished by their energy and particle fluxes. The weakly turbulent process can be 
characterized by fluxes in wavenumber space, while additional fluxes in amplitude 
space emerge in the intermittent process. The particle flux is directed from low am- 
plitudes towards high amplitudes, and the energy flows in the opposite direction. 
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1 Introduction 



This paper investigates two different types of turbulence, firstly, the energy 
transport of weakly interacting dispersive waves with low amplitudes [1], sec- 
ondly, an intermittent type of turbulence [2] where localized high-amplitude 
structures appear from a low-amplitude background, and decay rapidly [3] 
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[4] [5] . Important cases of weak turbulence have been found in the dynamics of 
Langmuir waves in plasmas [6] , capillary [7] and gravity waves [8] in hydrody- 
namics, in magnetohydrodynamics [9], and in Bose-Einstein condensates [10]. 
Collapses in Langmuir- waves [11] and self- focusing effects in optics [4], [12] are 
examples for the intermittent kind of behavior. 

Turbulent energy flows emerge when such dynamical systems are driven on a 
long spatial scale, and damped on a short viscous scale [1]. Models of turbulent 
flows assume that the dynamics is almost Hamiltonian in the inertial range 
between the damping range and the driving range, and energy is transported 
from sources to sinks through this window of transparency with no damping 
or driving. Symmetries besides the time-invariance are related to conserved 
quantities in addition to the Hamiltonian. Examples for such quantities are 
the light intensity in nonlinear optics [12] or the average momenta in hydro- 
dynamics [13]. When supplied and dissipated by driving and damping, there 
are fluxes also of these quantities. Fluxes that are directed from a source both 
toward long and toward short length scales can be modeled with an additional 
sink at small wavenumbers. 

The fluxes of the conserved quantities determine the statistically stationary 
nonequilibrium of the Hamiltonian dynamics in the inertial range. The theory 
of weak wave turbulence computes the nonequilibrium state of weakly inter- 
acting dispersive waves, depending on boundary conditions for the fluxes. A 
closed kinetic equation for the two-point functions can be derived if, with the 
existence of a small parameter, the nonlinear coupling is much smaller than 
the linear dispersive effects, so that higher order correlations can be reduced 
to powers of the two-point function. The kinetic equations can be solved an- 
alytically. This yields the Kolmogorov-Zakharov spectrum and formulas for 
the fluxes [1]. The theory predicts that the system will have close to Gaus- 
sian statistics even after the weak nonlinearity is relevant. The intermittent 
occurrence of high-amplitude structures [14] violates the assumptions of weak 
turbulence, as these structures indicate strongly coherent waves. Consequently, 
the closure of the kinetic equation is not possible, and there is a strongly non- 
linear interaction. 

In recent years massive numerical experiments have been carried out [15], 
[16], [17], [18] to check the assumptions of the Gaussian closure of the kinetic 
equations. This was partly thanks to the development of dynamical system by 
Majda, McLaughlin, and Tabak ('MMT' equation) [15], [16], which we use in 
its most simple focusing version 



.dcj) d 
dt dx 



1,2 <p-<p\<p\ 2 + F + v. (1) 



is a complex field that depends on one spatial variable x. The fractional 
derivative Id/dx] 1 ^ 2 ^ stands for the Fourier-transform of the term ^J\k\ak with 
the mode au = (l/£) / (f>(x) exp (ikx)dx, or equivalently, for the inverse of 
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a fractional integral [19]. The concave dispersion u(k) = y\k\ yields more 
four wave resonances than a convex dispersion law (such as the dispersion 
u(k) ~ k 2 of the nonlinear Schrodinger equation) does. This leads to an effi- 
cient energy transport even in only one dimension, and allows numerical sim- 
ulations in a wider range in the wavenumber space compared to simulations 
in two or three spatial dimensions. The general MMT equation has fractional 
derivatives of any order both in the linear and in the nonlinear part, but the 
semi derivative and the simple nonlinearity of (1) is by the signs of the fluxes 
[17] suitable for the study of wave turbulence. 

V is a damping term that affects very short waves and very long waves only. 
The driving force T affects waves on an intermediate length scale. The length 
scale of the driving force and the length scale of the short-wave damping are 
widely separated by the window of transparency where no damping or driving 
occurs. 

The MMT equation for zero driving JF and damping T> derives from the Hamil- 
tonian (or "energy") E = E 2 + E 4 , with the quadratic coupling part E 2 = 
J Wd/dx^^^dx = £ Vk\a k \ 2 , and the nonlinear part E 4 = - J \<j)\ A /2dx. The 
wave action or "particle number" n = J \<fr\ 2 dx = J2 |&fc| 2 is a second conserved 
quantity which is related to the system's phase symmetry. 
Numerical simulations of (1) found the Kolmogorov-Zakharov spectrum as 
predicted by weak turbulence theory, but also concluded that most of the 
energy transport is related to strongly nonlinear collapses [17], which coex- 
ist with weakly interacting low-amplitude waves. In this case, high-amplitude 
peaks emerge in a very similar way as in the corresponding Hamiltonian sys- 
tems where this process maximizes the entropy [20] . Pure wave turbulence for 
the MMT equation has been found for freely decaying turbulence [16], but not 
in the stationary nonequilibrium. A steeper spectrum than predicted by wave 
turbulence has been obtained for the defocusing MMT equation (where the 
sign of the nonlinearity is opposite to the focusing equation (1)) and in some 
cases for the focusing equation [16]. This lead to the proposal of an alternative 
closure from which the steeper spectrum can be derived [15]. 
The aim of this paper is to study both the weakly turbulent state and the 
collapse dynamics by measuring the fluxes of energy and particles. In fact, as 
we show here, flux measurements are a more reliable indicator of the transport 
dynamics than the energy spectrum. We compute energy and particle fluxes 
by measuring the input and output of these quantities that is caused by the 
driving force and the damping in our numerical experiments. 
We give evidence that a pure weakly turbulent process emerges if the driving 
forces are sufficiently weak. The fluxes in wavenumber space correspond to the 
predictions of weak turbulence theory (section 3). 

For stronger external driving we find a transition to a process where weak 
turbulence and collapse dynamics coexist. Despite the fact that the forma- 
tion of high-amplitude structures is not in agreement with the assumptions 
of weak turbulence our simulations of (1) find the weak turbulence spectrum. 
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The collapse dynamics is characterized by fluxes of energy and particles both 
in wavenumber and in amplitude space. These fluxes in amplitude space con- 
tain information on the coherent process (section 4). 

We present numerical results for one simulation that is representative for weak 
driving and damping (Figs. 2,3,4) and one for strong driving and damping 
(Figs. 2,6,7,8). Figs. 9,10 show data from another simulation with a medium 
power driving force. These two figures show time- dependent quantities, all 
other figures show time-averaged quantities. Figs. 5,11 show results of a series 
of simulations with various driving forces. 



2 Numerical method 

We integrate equation (1) numerically using a pseudospectral method with an 
integrating factor similar to the method described in [15]. <f>(x) is defined on 
the interval [0, 4096] with periodic boundary conditions. For its integration we 
use L = 4096 modes, with wavenumbers k in the Brillouin zone ] — 7r, 7r]. The 
linear part can be eliminated analytically using an integrating factor, and the 
resulting non-stiff equation is integrated with a multistep (Adams) method. 
We use driving and damping forces that act as time-periodic 5-kicks, while 
the system is governed by the Hamiltonian MMT equation for one time unit 
between the subsequent kicks. The external driving force T is applied to the 
modes in the interval 20 < \k\L/(2n) < 40. The driving 5-kicks increase the 
amplitude as a k — > (1 + \)a k with A = 0.1. This driving force is only applied 
to modes with an amplitude below a threshold |afc| 2 < co/k, so that the total 
energy input reduces when more modes exceed this threshold. 
The damping V is applied both at low wavenumbers (\k\L/(2n) < 20) and 
at high wavenumbers n/2 < \k\, while modes in the window of transparency 
40 < \k\L/(27i) < L/4 are neither damped and nor driven. The broad high-/c 
damping minimizes the aliasing error of the numerical computation. Ideally, 
the low-/c damping and driving range should also be chosen larger, but this 
would narrow the window of transparency. The damping 5-kicks decrease the 
amplitude as a k — > (1 — X)a k with A = 0.1. This damping is applied only 
to modes in the high-A; window when the amplitude exceeds the threshold 
|afc| 2 > c + /\k\, and to low- A; modes above the level |afc| 2 > c_. This type of 
damping becomes more effective when more modes exceed this threshold. 
The purpose of this type of damping and driving is purely practical: Damping 
and driving by time-periodic 5-kicks has the advantage that the input and out- 
put of energy and particles can be measured directly, and are distinguishable 
from the Hamiltonian dynamics of (1). The effect of the kicks is quite small in 
all our simulations, so that there is little difference between snapshots before 
and after the kicks. The damping is restricted to the two windows in £>space 
in order to keep the window of transparency free of dissipation. The threshold 
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Fig. 1. The driving force acts on modes with 20 < \k\L/2ir < 40 when [a^) 2 is 
below a threshold co/k. The driving (5-kicks increase these modes as a& — ► (1 + A)a^ 
with A = 0.1. The damping decreases the modes of the interval \k\L/2ir < 20 as 
flfc - ► (1 — A)afc when |afc| 2 above a constant threshold c_, as well as the modes in 
the interval 7r/2 < \k\ < tt when |afc| 2 is above a threshold c + //c. 

parameters c_, Cq and c + allow us to control the amplitudes in the driving and 
the damping regions. In particular, the driving force can be increased by in- 
creasing the threshold c . This type of forcing is similar to coupling the system 
to heat baths with different temperatures in the driving and damping regions. 
It appears that this leads to a relatively quick relaxation to a statistically 
stationary state. In this state, it is equivalent to other types of forcing that 
cause the same fluxes. In our simulations we choose c_ = Co/60 and Co slightly 
larger than c+. This choice has turned out to yield Kolmogorov-Zakharov 
spectra with relatively small deviations at the boundaries of the window of 
transparency. 

For comparison, we have performed simulations with time-continuous driv- 
ing and hyperviscous damping, which show no significant differences to the 
5-driven case. Our results are also in agreement with those reported in the lit- 
erature. In particular, we have found WT spectra [17] for the focusing MMT 
equation (1), while we find the steeper MMT spectrum [15] for its defocus- 
ing version. The results in this paper are robust under various modifications 
of the scheme. We have performed simulations with various resolutions from 
1024 modes to 8192 modes, and various types of damping and driving. Other 
modifications include the ranges in /c-space where damping and driving is ap- 
plied, the time intervals between the kicks and the strength of the kicks. We 
have tested of the integration scheme for undamped and undriven equations 
by simulating the well-known solitons of the nonlinear Schrodinger equation, 
and by checking the conservation of energy and particles by the MMT equa- 
tion. 

In all simulations the system is first allowed to relax until it reaches an almost 
stationary nonequilibrium state, which takes 10 4 — 10 5 time units depending 
on the particle input. Time-averaged quantities such as power spectra and 
energy fluxes are computed for the nonequilibrium state in a second round 
of integration. In this state, the input of particles and energy matches the 
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output, and the energies E% and E± and particle number n are constant with 
only small fluctuations. The integration time is long enough to average away 
the effect of these random fluctuations. 



3 The weakly turbulent state 

3. 1 Power spectrum and amplitude statistics 

The dynamics of (1) depends critically on the strength of the driving force and 
the damping, and we study the dynamics in two representative simulations, 
one for weak driving, and one for strong driving. Fig. 2 (a) shows the time- 
averaged wave action rik =< \ak\ 2 + \a^k\ 2 > over the wavenumber. Both for 
the weakly and the strongly driven system, we obtain a Kolmogorov-Zakharov 
spectrum rik = c/(kL/2ir) in the window of transparency. The coefficient c 
measures the particle number for waves in the window of transparency, which 
is c = 2.4 for the weakly driven system and c = 27 for the strongly driven 
system. In the regions in /c-space where the driving (damping) is applied, rik is 
higher (lower) than the level of the spectrum rik = c/(kL/2n). In the weakly 
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Fig. 2. (a): The time-averaged wave-action rik =< |afc| 2 + | cx A; 1 2 > (n = J2 n k is the 

total particle number) of the MMT equation exhibits Kolmogorov-Zakharov spectra 
nk ~ k~ l in the window of transparency both for weak and for strong driving forces. 
The dotted lines are fits n& = c/(kL/2n) with c = 2.4 and c = 27 of the spectra 
in the window of transparency, (b): Probability density of the amplitude \4>\ for the 
weakly and the strongly driven system. The probability for high amplitudes decays 
as ~ exp(— \cf)\ 2 ) for the weakly driven system. The shoulder-shaped distribution 
for the strongly driven system exhibits a significant probability for high-amplitude 
structures (\4>\ « 1.5). Particle number and energy fluctuate about constant values 
in both simulations (weak driving: n w 18, E2 ~ 7.5, £4 f» —0.08, strong driving: 
n w 150, E 2 w 75, E A w -10). 
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driven system, the spectrum has also two humps at harmonics of the driving 
force at the low- A; end of the window of transparency. These effects could be 
reduced with a driving force that acts on a broader /c-band. 
In Fig. 2(b), p(|0|)<i|0| gives the probability that the amplitude at one location 
is between |0| and \<f>\ + d\(f>\. The distribution of the amplitudes is qualita- 
tively different for the two simulations. In the weakly driven system we find a 
distribution p ~ \<fi\exp(— \<fi\ 2 ) without any incident of an amplitude |0| > 0.3 
during 5 • 10 5 integration cycles. The probability density of the strongly driven 
system is again gathered at small |0|, but its shoulder-like shape shows that 
amplitudes of the order of one are obtained with a probability of about 1CT 3 . 

3.2 Input and output of particles and energy 

Fig. 3 shows the input and output rates of particles and energy by the driving 
and damping forces as functions of k for the simulation with weak driving 
of Fig.2. For clarity, low wavenumbers (a) and high wavenumbers (b) are 
plotted on different scales. Positive An k show the input of particles in the 
driving range (Fig. 3 (a)), while negative An k indicate particle losses both in 
the \ow-k (Fig. 3 (a)) and in the high- A; damping ranges (Fig. 3(b)). The total 
particle number is unchanged since input and output of particles match as 
J2 AnkAk = 0, and the particle number is also conserved by the Hamiltonian 
dynamics. 

A change An k of the particle density at the mode k also changes the density 
of the coupling energy Ei at this wavenumber by u{k)AnkAk. Consequently, 
the gains and losses of particles by the driving and damping translate to gains 
and losses of energy via the coefficient uo(k) = \fk. The energy gain in the 




(a) (b) 

Fig. 3. Particle- and energy-balance for the weakly driven system of Fig.2. The 
driving and damping forces change the particle contents of the mode k by An^AA; 
per time unit. Anj is the change of the particle density, Ak = 27r/4096 is the 
distance between adjacent modes. This particle input or output is connected to a 
change of the coupling energy by oj(k)AnkAk. (a): The low-A; damping range and 
driving range, (b): The window of transparency and the high- A; damping range. 
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driving range matches the energy loss in the two damping ranges, so that 
J^^ityArikAk = 0. The influence of driving and damping on the nonlinear 
energy E4 is very small. 

3. 3 Particle flux of the weakly driven system 

In a statistically stationary nonequilibrium state, the particle number \ak\ 2 + 
|a_fc| 2 at I A; I fluctuates about a time-independent average value rik (Fig. 2a). 
Particles that are injected in the driving range flow to the two damping re- 
gions, where particles are removed at the same rate. This defines two fluxes of 
particles, one in a direct cascade through the window of transparency (where 
no particles are injected or removed) to the high-A; damping, and one in an 
inverse cascade to the low-A; damping range. The particle flux q(k) through 
/c-space follows from the continuity equation for the particles 

Aq(k)/Ak - An k / At = (2) 

where AA; = 27r/4096 and At = 1 in our simulations (a corresponding par- 
tial differential equation A — > d applies to continously damped and driven, 
infinite systems). Arik is the particle input or output by the external driving 
and damping, which is minus the particle change that is caused by the flux 




Fig. 4. Particle- and energy-Muxes in cj-space for the weakly driven system of Fig. 2. 
(a): The particle flux Q is positive (i.e. particles move toward high uj) in the window 
of transparency and in the high-w (or high-A;) damping range, and it is strongly 
negative in the 1ow-cj damping range and in most of the driving range. Its gradient 
is negative in the damping ranges, where particles are removed from the system, 
and positive in the driving range, where particles are injected into the system. Q 
is constant in the window of transparency. It is zero at to = and at uj = ^/tt, 
as the particle input matches the particle output, (b): The energy flux P has the 
same qualitative characteristics as Q. Quj is the energy drift of the particle flux. 
K = P — Quj is the energy transfer from the inverse cascade at low uj (where energy 
is released) to the direct cascade where energy is absorbed by the energy drift. 
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q(k). This definition yields the minus sign in equation (2). The particle flux 
is obtained as q(k) = Z)| =0 An^Ak/At. As the system's particle number is 
constant, the particle flux vanishes at k = and at the short-wave boundary 
of the Brillouin-zone k = n. The flux can be defined in frequency space as 
Qiyj) = q{k{uj)) which has similar properties since uj{k) = \J\k\ is bijective 
in the interval [0, 7r]. Fig. 4 shows the particle flux Q(uj) for the simulation of 
Fig. 3. A negative sign indicates a flux toward small frequencies (or wavenum- 
bers), a positive sign shows a flux to high uj. The flux is strongly negative 
in the the inverse cascade at low uj. The sign changes at some point in the 
driving region, and in the window of transparency there is a small positive 
flux to high uj. The slope is positive in the driving range, and negative in 
the damping ranges. The flux of this direct cascade is /c-independent in the 
window of transparency, and decays in the damping region at k > it/ 2 where 
particles are absorbed. 

3.4 Energy flux in the weakly driven system 

Similarly to the particles, the energy that is injected at the driving scale flows 
to the two damping ranges. The energy flux p(k) = Z)jL u(k)An~ k Ak/At can 
be defined from the energy continuity equation 

Ap(k)/Ak - uj(k)An k /At = (3) 

where the energy sources and sinks are located at the same wavenumbers as 
the particles sources. The constant energy of the nonequilibrium state requires 
that there is no flux through the boundaries p(0) = p(ir) = 0. Fig.4(b) shows 
the corresponding flux in u;-space P{uj) = p{k{yj)). 

As particles move toward higher values of uj in the direct cascade, their load 
with coupling energy increases proportionally to uj. The energy transport of 
the particle flux is measured by the quantity Q(uj)uj, which we call the energy 
drift. The energy drift depends on uj even in regions that are remote from 
energy sources and sinks: In the window of transparency, it increases with 
uj while the energy flux P(uj) and the particle flux Qiyj) are constant. The 
inverse cascade is linked to a negative energy drift that approaches zero for 
uj -> 0. 

K = P — Quj amounts for the discrepancy of the energy flux and the energy 
drift. It measures the part of the energy flux P that is not related to the 
particle flux. We call K the energy transfer term, because it describes the 
transfer of energy from the inverse cascade to the direct cascade: Particles 
of the inverse cascade release energy as they move towards lower uj, which 
feeds the energy transfer process. The energy transfer K shifts this energy 
toward higher frequencies by four-wave interactions. This energy is absorbed 



9 



by the energy drift in the direct cascade, where the particles demand energy 
in order to move to higher oo. In other words, the inverse cascade drives the 
direct cascade via the energy transfer. K is positive for all wavenumbers in 
this simulation, and vanishes at the boundaries. 

Fig. 4 and Fig. 2 give evidence for a pure weakly turbulent flow without strongly 
nonlinear processes. The amplitudes are always small, and the coupling part 
of the Hamiltonian carries almost all energy. The nonlinearity acts only as a 
weak coupling force of linear waves, without directly interfering in the energy 
input or output. As suggested in [17], the main ingredient to obtain such a 
solution numerically is a sufficient number of modes in the low- A; damping 
range and in the driving range. This avoids a bottleneck in the particle flow 
of the inverse cascade due to the lack of resonant modes. 
The Zakharov solution n k ~ 1/k to the kinetic equation involves a finite 
energy flux P > and a zero particle flux Q = in the direct cascade. This 
may be regarded as the limit of the fluxes of Fig. 4 for an infinite extension of 
the window of transparency, with a high-frequency damping at uo — > oo. The 
energy transfer K is again fed by the release of energy in the inverse cascade, 
but it is qualitatively different in the direct cascade. The Zakharov solution 
has a zero energy drift in the direct cascade, so that the energy transfer equals 
the energy flux K — P, and is ^-independent in the window of transparency. 

3. 5 Transition of the particle flux 

We now determine how the strength of the particle flux changes with the 
particle number. We run a number of simulations with various driving powers 
and determine the coefficient c of the spectrum n& = c/(kL/27r) for each of 
these simulations, c = 2.4 of Figs. 3,4 corresponds to the weakest driving force. 
Fig. 5 shows the strength of the particle fluxes Q of the direct cascade and 
the inverse cascade for various driving forces. Q + is the particle flux in the 
window of transparency, where the particle flux of the direct cascade reaches its 
maximum. \Q-\ is the modulus of the particle flux at the boundary between 
the low- A; damping range and the driving range (see Figs. 3 (a), 4(a)). Both 
fluxes increase with c, and the lines in Fig. 5 (a) are tentative interpolations 
Q ~ c 3 - 45 . For comparison, the wave turbulence closure gives Q ~ c 3 for four- 
wave interactions. 

The input of particles and of energy matches the output only for a certain ratio 
[1] of the strengths of the fluxes Q + and Q_. We assume for simplicity that the 
driving force acts only at one effective frequency cu , and that the high- and 
low-frequency damping act at u + and uj_ respectively. The particle output 
rates at these sinks are Q + and \Q-\, with corresponding energy output rates 
Q+uj + and \Q-\cu-. Input of particles Q + + \Q-\ and of energy (Q + + |<3-|Vo 
and their output match if Q+/\Q-\ = (ou — c<j_)/ (u + — u ) m 0.05, where we 
have set oo + xs 1.5, lu_ 0.15, uj ks 0.21. Fig. 5(b) shows this ratio of the 
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Fig. 5. (a): Maximum particle flux of the direct cascade Q + and of the inverse 
cascade \Q-\ versus the coefficient c of the wave action n& = c/(kL/2ir). The lines 
give a fit with ~ c 3 45 . (b): Ratio of direct and inverse cascade versus the wave action. 
For low particle densities, the ratio agrees to the prediction of weak turbulence 
theory Q+/IQ-I = ( w o — UJ -)/( u} + ~ ^o) ~ 0.05, where uj + rs 1.5 is the average 
frequency of high-/c damping and u;_ 0.15 of low-A; damping, while ujq 0.21 is 
the frequency where the driving occurs. 

two fluxes. For weak driving forces (corresponding to c < 10), Q + and \Q_\ 
both increase at the same rate with c, so that we find the expected fixed ratio 
of fluxes Q + /\Q-\ ~ 0.05 (line in Fig. 5(b)). Above c ~ 10 however, the flux 
\Q-\ increases at a lower rate with c, and the ratio Q + /\Q-\ increases. This 
is obviously not compatible with the assumption of a weakly turbulent state. 



4 Collapses and fluxes in amplitude space 



4-1 Particle flux and energy flux for the strongly driven system 



Above the transition point of Fig. 5(b) we observe intermittent high-amplitude 
structures. These structures are visible in the amplitude statistics of Fig. 2(b), 
but not in the power spectrum Fig. 2 (a) since only few of the particles con- 
tribute to collapses, and their Phillips spectrum decays more rapidly with k 
than the WT spectrum of the low-amplitude background [18]. We focus on 
the fluxes of particles and energy related to the strongly nonlinear events, and 
mention their dynamical properties only briefly. 

Figs. 6-8 show results for the strongly driven system of Fig.2 with c = 27. 
Fig. 6 shows the fluxes of particles and energy in wavenumber space. Again as 
in Fig. 3(a), the particle flux in Fig.6(a) is negative at small oj and positive 
above the driving range, and it reaches zero both for high and for low uo. Com- 
pared to the weakly driven system of Fig. 3, the particle fluxes are stronger by 
three orders of magnitude, and a higher share of particles (one fourth com- 
pared to one twentieth) moves toward high u. The relative strength of the 
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Fig. 6. Particle- and energy-fluxes of the strongly driven system of Fig. 2. (a): The 
particle flux Q has the same shape as for the weakly driven system (Fig.3), but 
the positive flux towards high uj has a higher relative strength. Again there is no 
net gain or loss of particles, (b): Unlike for the weakly driven system (Fig. 3b), the 
energy flux P and the energy transfer K become negative for high to. 

direct cascade leads to a significant change of the energy fluxes in Fig. 6(b): 
The dissipation at high ui causes to a decay of the energy flux P below zero. 
The dissipation of coupling energy exceeds its input by AE 2 within each cycle 
of driving and damping. Contrary to this apparent loss of coupling energy, E 2 
is constant with only small fluctuations, and consequently there must be an 
additional source of coupling energy to compensate this loss. 
This source also plays a role for the energy transfer K. K is positive in the 
lower part of the window of transparency, which indicates an energy transfer 
from the inverse cascade that is feeding the energy drift in the direct cascade. 
However, the sign of K changes within the window of transparency, and K < 
suggests that energy is transfered from a source at high ui to lower frequencies, 
where it feeds the energy drift Qui. The energy drift is therefore now driven 
both by the energy transfer from the inverse cascade, and from a transfer of 
energy from some other source in the high-cj dissipation range. We will now 
study this source of energy. 



4-2 Nonlinear energy 

Fig. 7 shows the average input and output of particles and of nonlinear energy 
£4 as a function of the amplitude |0|. For a formal definition, we define the 
ceiling function 

{v for v < w 
w for v > w 
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for real v, w. The square norm of the ceiling function of the amplitude \4>(x) | is 
an amplitude-dependent particle number n(|0|) = / c 2 (|0(x)|, \(f)\)dx. This par- 
ticle number is obtained when the amplitude is set to equal |0| at all locations 
where it exceeds |0|. The derivative n\M = dn(\(f>\) / d\(f>\ gives the nonnegative 
density of particles as a function of the amplitude |0|, and n = Jo° n\Md\4>\ * s 
the total particle number. An\M as shown in Fig. 7(a), (b) is the time-averaged 
difference of jiui before and after the driving and damping kick. In other words, 
jiui is the particle density in amplitude space that is injected or removed by 
the driving and damping. There is a net particle input at low amplitudes and 
a net particle output at high amplitudes. The area under the curve where par- 
ticles are gained equals the area under the curve where particles are removed, 
as there is no net gain or loss of particles. 

Analogously to the particle density, a density of nonlinear energy in ampli- 
tude space Emm = dE^Qcf)]) / d(\(f>\) derives from the amplitude-dependent en- 
ergy E 4 (\(j)\) = -/c 4 (|0(a;)|, \(j)\)/2dx. AE^\ of Fig.7(a),(c) is the average 
difference of E^ before and after the driving and damping. This shows a 
net gain of nonlinear energy at high \<f>\ where particles are dissipated, and a 
small net energy loss at low \(f>\ where particles are injected. As the change 
of the amplitude A|0| <C |0| is small compared to the amplitude, the change 



2 driving 




Fig. 7. Balance of particles n and nonlinear energy E4 as a function of the amplitude 
for the strongly driven system of Fig. 2. (a): Damping and driving lead to a net 
particle gain at those locations in real space where the amplitude is small (|^»| ~ 0.3) 
and to a net particle loss at higher amplitudes (|0| w 1.4). The net losses match 
the net gains. The losses at high \4>\ lead to a gain of nonlinear energy — |0| 2 Ani^i. 

(b) : Driving and damping inject and remove particles mainly at low \<p\. The net 
particle change (see also (a)) is the sum of the driving and the damping process. 

(c) : The net change of the energy E4 (see also (a)) as the sum of the driving and 
the damping process. 
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of the energy density is connected to the change of the particle density by 
AE^m « — |0| 2 Anj^|. Sources and sinks of the particles translate to sinks and 
sources of nonlinear energy via the factor — 10| 2 . This is somewhat similar to 
the correspondence between sources of particles and coupling energy in /c-space 
of Fig.3, where the coefficient is uj(k). The particle sink at high \(f)\ becomes a 
substantial source of energy E4, as the factor — 10| 2 is quite significant there. 
At the same time, the particle source at low \<f>\ is only a small energy sink. 
Figs. 7 (b) and (c) distinguish between the effects of driving, damping, and 
the net effect of the superposition of both (as in (a)). Fig. 7 (b) shows that 
that most particles are both injected and dissipated at locations where the 
amplitude is small. The sum of both curves shows again the net particle input 
at |0| ~ 0.3 and the output at \(f>\ ~ 1.4. Fig.7(c) shows the corresponding 
changes of nonlinear energy, where the gain at high amplitudes is the most 
important effect. In contrast to this, no relevant change of the nonlinear en- 
ergy is observed in the simulation with weak driving of Fig.3. In this case, the 
particles are both injected and removed at low amplitudes, and there is no 
particle dissipation and no energy gain at high \<f>\. 

4-3 Fluxes in amplitude space and energy from collapses 

The particle source at low \cf>\ and the sink at high \<j)\ induce a particle flux 
A(|0|) through amplitude space. Fig.8(a) shows this particle flux, which is 
obtained from dA(\(j>\)/d\<j>\ — An^/At = by integration of An^| over \<f>\. 
Similarly to the particle fluxes Q in cu-space (Figs.3(a),4(a)), it is zero both 
for low and for high |0|, as the total particle number is constant. A(|0|) is 
positive, which shows that the particle flux is directed from low amplitudes 
to high amplitudes. Unlike the fluxes in k or w-space, the sources and sinks 
overlap, so that only a part of the total particle input contributes to the flux, 
while most particles are again removed by the damping at the same |0|. 
A nonlinear energy "flux" M(|0|) is obtained by integrating <9M(|0|)/<9|0| — 
A-Ejui/Ai = 0. This quantity is negative at small amplitudes, since the par- 
ticle sources at low |0| are sinks of £4. At higher |0| it increases with the 
amplitude, since the particle sinks in this domain are powerful sources of £4. 
It reaches a saturation value A£ 4 at |0| ~ 1.8, which is the net gain of non- 
linear energy by the driving and the damping. 

In a similar way, a "flux" II(|0|) of coupling energy in amplitude space can be 
defined. The coupling energy of the locations x where the amplitude <ft(x) is be- 
low the threshold \ij>(x)\ < |0| is E 2 (\(f)\) = f u(\(j>\ - \^(x)\)\\d / dx\ l l^(x)\ 2 dx, 
where u(l > 0) = \,u{l < 0) = is the unit step function. The corresponding 
energy density in amplitude space is £2Ui — dE2{\<p\) / d\<p\. We measure the 
time-averaged energy input or output A-E^ui of the density of the coupling 
energy in amplitude space. II(|0|) is obtained by integrating Ai? 2 |0| over \<p\ 
(Fig. 8(d)). II(|0|) is zero at |0| = 0, and decreases to a constant value AE 2 < 
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Fig. 8. Particle and energy fluxes in amplitude space, (a) Particle flux A(|(/>|) as 
obtained by integration of the net particle balance of Fig. 7. Particles move towards 
high \(p\ as coherent structures grow. A(|^>|) is zero both for \(f>\ = and for high \(f>\, 
as the system neither gains nor loses particles, (b) The "flux" M(|^|) of nonlinear 
energy, (c) The total flux of energy (nonlinear energy M(|(/>|) plus coupling energy 
I1(|0|)) is negative, i.e. energy flows from high amplitudes to low amplitudes, (d) 
The "flux" of coupling energy n(|0|) is negative. The loss of coupling energy AE 2 
in the flux n(|0|) is identical with the loss of coupling energy in the flux P in Fig. 4. 
The loss of coupling energy AE 2 matches the gain AE4 of nonlinear energy. 

at \4>\ ~ 1.8. n(|0|) expresses in amplitude space the same loss of coupling en- 
ergy as Fig. 6(b) does for the flux P(w) in frequency space. This loss of coupling 
energy AE 2 matches the gain of nonlinear energy AE^ of Fig.8(b). 
These gains and losses are caused by the external driving and damping, but 
the Hamiltonian dynamics reverses these changes. £4 and E 2 each are not 
conserved under the Hamiltonian dynamics, which actually transforms non- 
linear energy into coupling energy at the same rate as damping and driving 
feed in nonlinear energy, and dissipate coupling energy. As a consequence, £4 
and E 2 each have no gains or losses after an average cycle of driving, damping, 
and Hamiltonian evolution over one time unit. Nonlinear energy is produced 
mainly by the dissipation of particles at locations where the amplitude is high. 
The dissipation of coupling energy is also relevant for low-amplitude fluctua- 
tions. The total energy flux M(|0|)+n(|0|) (Fig.8(c)) is therefore directed from 
high amplitudes to low amplitudes. In wavenumber space, this is related to 
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the negative energy transfer K(u) at high oo (Fig. 6(b)). K is fed by the source 
of nonlinear energy, which is due to the dissipation at high frequencies. This 
energy is transformed into coupling energy in the window of transparency. 
M(|0|), n(|0|) and P{oo) are not genuine fluxes, because E 2 and E 4 each are not 
conserved under the Hamiltonian dynamics. This causes the false impression 
of Figs. 6(b), 8(b), 8(d) of permanent losses of E 2 and gains of E4. These quanti- 
ties do not account for the k- nonlocal energy transformation of £4 into E 2 . In 
contrast, A(|0|) and Q{uj) as well asM(|0|)+Il(|0|) are fluxes of the conserved 
quantities n and E respectively. In the weakly driven system (Fig.4(b)), P(u>) 
may also be considered as a flux, because E 2 is virtually constant under the 
Hamiltonian dynamics. 

4-4 The energy cycle of collapses 

We now investigate the Hamiltonian process that transforms energy E4 into E 2 
by studying the time-evolution, instead of averaged quantities. Figs. 9, 10 show 
a simulation with a medium-power driving force with c = 13.6. Fig. 9(a) shows 
the maximum square amplitude max(\(f)\ 2 (x)) that is reached at any location 
x in the system. The system has almost everywhere low amplitudes |0| 2 < 0.3, 
and high amplitude structures |0| 2 > 1 emerge only intermittently in small 
areas in space. The medium-power driving force allows us to distinguish dif- 
ferent amplitude blow-ups in time, while a stronger driving force would lead 
to overlapping high-amplitude signals from simultaneous blow-ups at different 
locations. The fluxes of the the cases c = 13.6 and c = 27 (Figs. 6-8) have 
similar qualitative features. The weak driving force with c = 2.4 of Fig. 4 leads 
to no such high-amplitude events, so that the maximum amplitude is of the 
order \<p\ 2 ~ 0.01. 

It is not completely clear how these structures emerge from the random low- 
amplitude background. The nonlinear events are very similar to collapses in 
Hamiltonian partial differential equations such as the undamped MMT equa- 
tion or the nonlinear Schrodinger equation. These can be described as het- 
eroclinic orbits that connect the phase-unstable homogeneous mode with a 
localized structure with an infinite amplitude. A collapse contains a finite 
amount of particles with a zero energy, since its nonlinear energy is minus 
the coupling energy. Collapses may be regarded as a coherent and /c-nonlocal 
mechanism that transports particles towards high wavenumbers [4], [17], [18]. 
Dissipation of short waves prevents the singularity, and leads to a decay of the 
peak. 

In the focusing nonlinear Schrodinger equation, all plane waves are phase- 
unstable under slow modulations, but this mechanism is different in the MMT 
equation. A perturbation 5(f) = exp(ikx)(a(t) exp(iqx) + b(t) exp(-iqx)) of a 
plane wave = y/nexp(i(kx — (u — n)t) has an eigenvalue with a positive real 
part if {2n — v)v > with v = (uj{k + q) +u{k — q) — 2u(k))/2, which requires 
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Fig. 9. Simulation with a medium-power driving force (c = 13.6). (a): Maximum of 
the squared amplitude \(f)\ 2 that is reached anywhere in the system as a function of 
time, (b) E2, E4 and E2 + E4 as functions of time, (c) Cumulated input of coupling 
energy by the driving force J2o AE^f™™ 9 ^ , and cumulated output of coupling energy 
by the damping J2o AE!f aTnpm9 ' > , (d) Cumulated change of coupling energy by the 
damping and the driving force Ya) AEif"^ 9 ^ + ^j^ am ~P m ^ ( sum Q f the curves of 
(c)), and cumulated change of nonlinear energy by the damping and the driving 
force £0 AE[ drivin9) + AEf amping) . 



< v < In. For the focusing nonlinear Schrodinger equation, this is always 
fulfilled for small q since u = k 2 is convex and v = q 2 > 0. For the MMT 
equation, a necessary condition for an instability is that k + q and k — q have 
different signs, i.e. |g| > \k\. The homogeneous mode k = is unstable under 
slow modulations, since the dispersion uj = ^J\k~\ is convex at this point. For 
waves with k 7^ 0, the dispersion is concave, and an instability occurs only 
by interaction with remote modes, i.e. modes on both the left and the right 
branch of uj. 

A possible starting point for collapses might be the built-up of a condensate 
at k = by the inverse cascade [4]. The phase instability of the condensate 
might lead to a collapse as soon as enough particles have gathered at k — 0. 
However, in our simulations the low-A; damping is so efficient that the homo- 
geneous mode has a very low level. Collapses in our simulations occur while 
the homogeneous mode contains only 11^=0 ~ 0.1 particles, i.e. the number of 
particles in the homogeneous mode of the whole system is below the particle 
number of a single collapse. In addition, we find no sign of a build-up of a 
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Fig. 10. (a): E 2 and E 4 as functions of time for the simulation of Fig.9. Energy 
is transferee! from £4 to E 2 during the growth of a peak. This energy is partly 
dissipated and partly transfered back to £4 when the peak decays. The damping 
increases the energy £4. (b): Energy gain / \^^ lam P m 9) anc [ ener gy j ss AE 2 dampmg ^ 
by the damping. 



condensate, as nk=o appears to be independent of the occurrence of collapses. 
Collapses occur in the same way for even more efficient low- A; damping, where 
each damping step removes all particles with small wavenumbers. The collapse 
dynamics appears to depend crucially on the maximum power in the driving 
range, or equivalently, on the total particle number. Instabilities of waves in 
the driving range are strongly nonlocal in wavenumber space and cause no 
slow phase modulations. 

Fig. 9(b) shows the energies E 2 , £4, and E = E 2 + E 4 of the total system 
during the simulation. At each of the high-amplitude events, £4 temporarily 
decreases, while E 2 increases by the same amount, and both energies approach 
their previous values subsequently. There is no net change of these energies 
on average. 

Fig.9 (c) shows the input and output of coupling energy as functions of time. 



(driving) 



is the input of coupling energy through the driving force cumu- 



lated over a period of time [0, t], and J^o AE 2 dampmg ^ is the corresponding out- 
put through the damping. Coupling energy is dissipated at an approximately 
constant rate between the collapses, which reflects a weakly turbulent trans- 
port of energy into the dissipation range between the collapses. The increased 
dissipation at times when high-amplitude structures emerge leads to step-like 



decays of £0 AE ( 2 



(damping) 



The steps lead to an energy output that is bigger 



than the energy input from the driving, i.e. £ A .Eg 



(damping) 



+ EoA£ 



(driving) 



< 

0. This net output of coupling energy is shown in Fig.9(d) together with the net 
input of nonlinear energy £0 A^f ampin9) + £* AE {d nving) > 0. The time aver- 
age of AE {dampmg) + AE {driving) is just the energy-flux P{u) = II(|0|) = AE 2 
at high 00 or high \(f>\. The average of [^j^amptng) ^^(dnvmg) ^ ^ e ener gy_ 
flux M(|0|) = AE 4 at high |0|. As there is no net change of the total energy, 
the gain of E 4 matches the loss of E 2 . 

Fig. 10 shows two of the collapses of the simulation of Fig.9. When a peak 
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Fig. 11. (a) Net loss of coupling energy AE 2 and gain of nonlinear energy AE4 
versus the coefficient c of the power law = c/(kL/2ir). AE 2 is the value of the 
energy flux P at high lo or II at high \<fi\, AE4 is the value of the flux of nonlinear 
energy M at high \<fi\. (b) Peak heights together with an upper and a lower estimate. 
The simulations are the same as in Fig. 5. 

emerges, the total coupling energy E 2 increases due to the high gradient of 
the peak. Simultaneously, the nonlinear energy £4 decreases because of the 
locally high amplitude. The blow-up transfers energy from E 4 to E 2 . 
The peak is strongly affected by the damping so that it looses coupling en- 
ergy. Fig. 10(b) shows this dissipation, where the spikes of the dissipation rate 
^■g(damping) correS p OI10 i the steps of J2o AE 2 dampmg ^ of Fig. 9(c). As the am- 
plitude of the peak decreases by the action of the damping, the nonlinear en- 
ergy increases again. Therefore the damping provides both an input of energy 
E 4 and an output of energy E 2 (Fig. 10(b)). The picture of the energy cycle is 
that the damping feeds nonlinear energy AE 4 into the system, that is trans- 
formed to coupling energy AE 2 and subsequently dissipated. The burnouts of 
the peaks are incomplete [4], and some of the coupling energy is not dissipated, 
but returned to the nonlinear energy by a Hamiltonian process. Radiation of 
the decaying peak might also contribute to the low-amplitude waves. 
Fig. 11 (a) shows the net output rate coupling energy AE 2 and the net input 
rate of nonlinear energy AE4 that is caused by damping and driving. As in 
Fig. 5, c measures the particle density that is obtained in the window of trans- 
parency for different driving powers. Again, we can distinguish between pure 
weak turbulence and collapses. For driving forces below the threshold c ps 10, 
the cumulated dissipation of coupling energy (as in Fig. 9(c)) has no steps, 
as there are no collapses. Average input and output of coupling energy are 



matching as AE 2 



(damping) 



^ ^(driving) 



AE 2 0, and the input of nonlinear 
is negligible. Above the thresh- 



energy AEf amping) ps and AEf riving) 
old, output of coupling energy and input of nonlinear energy by damping and 
driving match closely as AE 2 = —AE 4 . The lines give a tentative fit ~ c 4 . 
A characteristic peak height can be denned as the amplitude where the non- 
linear energy change by the damping and driving is maximal, i.e. \cf>\ of Fig. 7 
with the maximum |A_E 4 |^||. Fig.ll(b) shows this peak height \<f>\+ versus the 
coefficient c of the power spectrum. Again, high amplitudes appear only for 
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strong driving forces. Below the threshold, | AEmm | is very small and negligible 
in the energy balance. To give an estimate the height of these structures, we 
assume that An particles carrying an energy u; An are injected at cu , and 
they are removed at ui + where they carry the coupling energy u + An. The net 
energy loss AE 2 = (ou — u + )An from this process can be balanced by a gain 
of nonlinear energy due to a particle flow towards high amplitudes where the 
particles are dissipated. If we assume for simplicity that the damping reduces 
the amplitude \(f>\ by A|0| <C |0| in the interval Ax in space, an amount of 
An = 2 1 0| A | 0| Ax particles are removed and the nonlinear energy is increased 
by AE A = |0| 2 An. With the energy balance AE A + AE 2 = 0, we obtain 
|0| = — ujq as a lower estimate for the height of the peaks. For an upper 
estimate, we assume that the damping decreases the amplitude by A|0| = |0| 
and removes all particles An = |0| 2 Ax from this interval, so that the gain of 
nonlinear energy is AE 4 = |0| 4 Ax/2 = |0| 2 An/2. For this case, the energy 

balance gives |0| = \j2(cu + — co ) as an upper estimate of the peak height. 
If the particles input by the driving force has an average u ~ 0.2 and the 
particle output occurs at u + ps 1.5, we expect a peak height between |0| ss 1.1 
and |0| ss 1.6, compared to the numerical result |0| ~ 1.4. 



5 Conclusions 

We have studied the fluxes of energy and particles in a onedimensional tur- 
bulent system, the focusing Maj da-McLaughlin- Tabak equation. These fluxes 
reveal details about the statistically stationary state that cannot be observed 
directly in the conserved quantities, such as the spectrum of particles n& 
(Fig. 2). Fluxes in frequency space (Figs. 4, 6) have sources and sinks that are 
determined by the external driving and damping. Fluxes in amplitude space 
(Fig. 8) emerge spontaneously when the total particle input rate exceeds a 
threshold (Fig. 5). 

We have found two types of nonequilibria with similar spectra (Fig.2(a)), but 
very different energy fluxes. The nonequilibrium state for sufficiently weak 
driving is very similar to the weakly turbulent Zakharov solution. The aver- 
age amplitude is small, with an exponentially decreasing probability for higher 
amplitudes (Fig. 2(b)). Nonlinearity only plays a role as a weak coupling be- 
tween the dispersive waves, and no high-amplitude structures emerge. Particles 
that are injected at the driving scale flow toward the sinks both at high and 
at low wavenumbers (Fig.4(a)). We interpret the well-known quantities [1] Qu 
and K (Fig. 4(b)) as energy drift and energy transfer. In particular, K trans- 
fers energy that is released in the inverse cascade to the direct cascade, where 
it is is absorbed. The fluxes increase slightly faster with the particle density 
(Fig. 5 (a)) than predicted by the theory of weak turbulence. 
For stronger driving forces, high amplitude structures emerge intermittently, 
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and the interaction becomes strongly nonlinear. The energy balance changes 
in two ways: Firstly, more coupling energy is dissipated by the damping than 
supplied by the driving force. This leads to a negative energy flux P and en- 
ergy transfer K for high frequencies (Fig. 6(b)). Secondly, damping and driving 
provide a net surplus of nonlinear energy £4, which is due to the dissipation of 
particles at high amplitudes (Figs.7,8(b)). This is related to fluxes of particles 
and energy in amplitude space (Fig.8). There is a net particle flux from low 
amplitudes, where the driving force injects particles, toward high amplitudes, 
where particles are removed by the dissipation (Fig. 8 (a)). The dissipation of 
particles at the tip of high-amplitude structures is the source of an energy flux 
from high amplitudes to low amplitudes, where coupling energy is dissipated 
(Fig.8(c)). 

While driving and damping cause a permanent net input of £4 and an equal 
net output of E 2 , both the coupling energy and the nonlinear energy are 
constant with only small fluctuations (Fig.9(b). This is possible because the 
Hamiltonian dynamics transforms nonlinear energy into coupling energy within 
the energy flux from high to low amplitudes (Fig. 8(c)). This transformation of 
nonlinear energy into coupling energy occurs during the blow-ups (Figs. 9, 10), 
where the coupling energy increases and the nonlinear energy decreases. Cou- 
pling energy can be radiated into the surrounding low-amplitude waves, when, 
for example, a peak emerges out of two colliding solitary waves [21]. A peak 
that decays after reaching its maximum height can also radiate low-amplitude 
waves. The generation of coupling energy at high u feeds the energy transfer 
K, that is directed toward low frequencies at high uj. The representation of 
energy fluxes in frequency or wavenumber space is somehow misleading since 
it neglects the /c-nonlocal transformation nonlinear energy into coupling en- 
ergy. 

The nonlinear energy £4 is related to high-amplitude coherent structures, 
while the coupling energy E 2 has contributions both from these peaks and 
from low-amplitude fluctuations. The advantage of studying these two quan- 
tities (and not peaks and fluctuations) is that their fluxes are measurable sta- 
tistical quantities that distinguish weakly turbulent and coherent processes. 
The fluxes in wavenumber or frequency space and in amplitude space may be 
seen as two projections of the complex highdimensional transport process. In 
wavenumber space, driving and damping is related to localized sources and 
sinks of the fluxes. The fluxes yield the boundary conditions for the nonequi- 
librium in the inertial range. This is appropriate for the pure weakly turbulent 
process, where the statistics is almost Gaussian and higher order correlations 
can be reduced to two-point functions. If the statistics is non-Gaussian, the ex- 
ternal forces inject energy to the nonlinear part of the Hamiltonian £4 (i.e. the 
four-point correlation) and remove energy from E 2 , the two-point correlation. 
The Hamiltonian dynamics is characterized by the transformation of energy 
£4 into E 2 , which is actually the transfer of a conserved quantity between 
correlations of different order. The representation in amplitude space covers 
this spontaneous coherent process, which appears as an energy flux from high 
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amplitudes to low amplitudes. The representations in wavenumber space and 
in amplitude space together cover the energetics of this process. The many 
open questions on the energy transfer from E4 to E 2 concern, in particular, 
the onset of this coherent behavior. Also, the effect of the coherent structures 
on the low-amplitude process by radiation is not well understood. 
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